Testing symmetries in effective models of higher derivative field theories 
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Higher derivative field theories with interactions raise serious doubts about their validity due to 
severe energy instabilities. In many cases the implementation of a direct perturbation treatment to 
excise the dangerous negative-energies from a higher derivative field theory may lead to violations of 
Lorentz and other symmetries. In this work we study a perturbative formulation for higher derivative 
field theories that allows the construction of a low-energy effective field theory being a genuine 
perturbations over the ordinary-derivative theory and having a positive-defined Hamiltonian. We 
show that some discrete symmetries are recovered in the low-energy effective theory when the 
perturbative method to reduce the negative-energy degrees of freedom from the higher derivative 
theory is applied. In particular, we focus on the higher derivative Maxwell-Chern-Simons model 
which is a Lorentz invariant and parity-odd theory in 2 + 1 dimensions. The parity violation arises 
in the effective action of QED3 as a quantum correction from the massive fermionic sector. We 
obtain the effective field theory which remains Lorentz invariant, but parity invariant to the order 
considered in the perturbative expansion. 

PACS numbers: 11.30.Cp, ll.10.Ef, ll.10.Lm 



I. INTRODUCTION 



The standard model has been extremely well tested at 
presently attainable energies, nevertheless it is in general 
regarded as an effective theory valid up to certain high 
energy scale at which small imprints of a more fundamen- 
tal theory can appear. This is assumed in part because 
the underlying theory could provide a solution to the 
ultraviolet divergences in quantum field theories, incor- 
porate gravity or explain other incomplete issues. With 
the aim to find some of the possible low-energy effects of 
the more fundamental theory, usually considered to be 
a quantum gravity theory or related to a unified theory, 
several extended field theory models with small modifica- 
tions have been constructed. These small modifications 
have been proposed mainly through the use of non-higher 
dimensional operators with some of them including the 
possibility of Lorentz and CPT symmetry violations [l[ . 
On the contrary, models containing higher dimensional 
operators have been less studied or attractive due to the 
many difficulties involved in their formulation even 
when they can be incorporated without any symmetry 
breaking. In spite of their drawbacks presently higher 
derivative field theories continue to be strongly motivated 
due to the insights they are believed to provide in the elu- 
cidation of the more fundamental theory. For example, 
the increase in the degrees of freedom of a self- interacting 
harmonic oscillator due to time derivatives of increasing 
order obtained in the calculation of the effective action 
[H, lij , may serve to evidence that we are approximating 
a more fundamental extended object which will become 
later formally defined as the wave functional in the quan- 
tum theory. 
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Higher derivatives were first considered in field theories 
to get rid of the infinities associated to point particles [f| ■ 
They can improve ultraviolet properties in quantum field 
theories @ and gravity , although for a contrary point 
of view, see B. Moreover, they have been obtained from 
string theory [9[ , non commutative theory [ic| , derivative 
expansions of the effective action 0, H| and have been 
used in electrodynamics [ll|, dark energy physics (l2l . 
Il3j , inflation [l4[ , as ultraviolet regulators [15[ , Lee- Wick 
models [llj], and some other contexts fl7j . 

Most of the problems with higher derivative field the- 
ories such as instability, causality violations, nonunitary 
evolution and the possible emergence of quantum states 
with negative-norm called ghosts states are intimately 
connected to the fact that the energy has no bottom. 
That is, in general the higher order Hamiltonian, the one 
producing temporal evolution, has an unusual part of the 
spectrum taking infinite negative values in addition to 
the infinite tower of positive energies. For the noninter- 
acting theory the negative-energy modes are decoupled 
from the positive ones, they evolve independently, which 
eventually introduces no harm into the classical or quan- 
tum field theory. However, when turning on the interac- 
tions, both excitations couple, giving origin to Feynman 
vertices in the quantum field theory with the possibil- 
ity of one particle decaying into the other through the 
process of creation and annihilation. Therefore, in com- 
pliance with energy conservation, an infinite number of 
positive and negative energy states are allowed and actu- 
ally favored to occur driving the system to an unlimited 
particle production. In consequence, the vacuum state in 
the Fock representation becomes rapidly unstable. One 
can try to overcome these problems by passing to an 
alternative realization where the negative-energy states 
are exchanged by negative-norm states, however spoiling 
the unitary evolution of the quantum theory. By using 
path integrals techniques, it is conceivable to have posi- 
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tive transition amplitudes that behave as perturbations 
over the ordinary ones, thereby taking very small depar- 
tures from unitarity at low energies [la]. Also, there are 
proposals where these problems can be circumvented by 
modifying the standard internal product in a VT sym- 
metric model fl9j . modifying the usual energy interpre- 
tation Ho|, or using BRST symmetry [2l|. Much of the 
instability problem and their related issues can be un- 
derstood by studying mechanical models for which we 
strongly recommend the Refs. 113 . |2S 
In many occasions [25|, [26|, l27l [28 



231, Hi]. 

291 ] one wants to 



consider higher derivative theories when they describe 
small deviations to an ordinary-derivative theory (ordi- 
nary in the sense that they contain no more than first 
time-derivatives in the Lagrangian). In that case an ad- 
ditional problem may arise since, no matter how weak 
the higher derivative operators are coupled, in general, 
the modified theory suffers an increase in the degrees of 
freedom [301 ] . which ultimately will depend on how much 
the theory is constrained. Therefore, a consistent per- 
turbative formulation requires us to retain the original 
number of degrees of freedom dictated by the ordinary- 
derivative theory and to reproduce only the dynamical 
sector of the higher derivative theory that is well defined 
when the higher derivative operators are taken to zero: 
the perturbative sector. In addition, in the presence of 
interactions these perturbative degrees of freedom should 
have a stable evolution. 

To illustrate some of the above problems, consider the 
Lagrangian 



c = —f^f^ + |e Q ^(n^ a )(a^ 7 ), 

which, by using the field redefinitions 
A» = -j=(A»+gF*»), 



A" 



V2 



F* 



can be rewritten as 



C = -,F^--A^- 2 



□ A", 



(1) 



(2) 



(3) 



(4) 



where F„ y = d^A v ~ d v A^ 
F*^ = ie A1Q/3 F a/ 3 the pseudovector dual field 



is the strength tensor and 
At this 



level we can emphasize that: 

(1) The Lagrangian consists in the sum of a massless 
and a massive term both associated to low-energy and 
high-energy degrees of freedom, respectively. The mass- 
less Lagrangian is the only one depending linearly on g 
so it must be parity odd. 

(2) The massive field contributes with negative energy 
to the system making the Hamiltonian unbounded from 
below. 

(3) Since the low-energy and high-energy fields are de- 
coupled from each other we still have unitarity and sta- 
bility in the theory [24]. 



(4) The mass of the high-energy field depends nonana- 
lytically on the parameter g, which goes to infinity when 
'/ 

(5) By turning on the interactions, for example, via 
a fourth power of the gauge field (A^ A^) 2 , the theory 
probably collapses giving all the stability problems we 
have mentioned. 

(6) In order to approach well the low-energy dynamics 
the nonperturbative massive mode has to be removed 
from the higher derivative theory which is not clear how 
to achieve by making field redefinitions. 

In this paper we study the perturbative method for 
higher derivatives theories developed in Refs. @, [3l| . 
The method allows to retain the low-energy sector of 
the higher derivative theory and to construct a positive- 
defined effective Hamiltonian order by order in the ex- 
pansion parameter g. The perturbative method has been 
thoroughly applied in [32l l33l. |34|. and, in particular, it 
has been implemented as a starting point to quantize 
Lorentz violating higher derivative field theories [34|. In 
this work we attempt to go a step further. From the 
observation that Lagrangian ([T]) describes a photon at 
low energies (taking the limit g — > 0), it sounds plausi- 
ble to recover electromagnetic symmetries by removing 
all the high-energy degrees of freedom from the higher 
derivative theory. Thus, the main goal in this paper is to 
test whether the symmetries of a higher derivative field 
theory can be modified by applying the perturbative for- 
mulation. It is well known that symmetries in an effective 
theory may depend on the scheme of approximation used 
to obtain the low-energy limit. For example, a direct an- 
alytical expansion in the Hamiltonian or implementing a 
order reduction treatment using the equations of motion 
on the higher derivative Lagrangian may lead to differ- 
ent results for Lorentz and other symmetries [HI, [36| . In 
order to implement and test the perturbative method we 
will consider the Lagrangian ([1]). 

The organization of the paper is as follows. In Sec. II 
we obtain the solutions to the equations of motion for 
the higher derivative Maxwell-Chern-Simons (MCS). We 
develop the Hamiltonian formulation with all the con- 
straints included and we exhibit the negative energies in 
the total energy spectrum produced by the massive field. 
In Sec. Ill we review both the Hamiltonian formulation 
for higher derivatives field theories and the perturbative 
method. Section IV is the main part of this work. There 
by implementing the perturbative method we derive the 
effective theory and we perform a complete study of its 
symmetries. In Sec. IV we give the conclusions and fi- 
nal comments. For completeness we provide the reduced 
Hamiltonian and the Dirac brackets in the appendix. 



II. THE HIGHER DERIVATIVE GAUGE 
THEORY 

In this section we will consider higher derivative cor- 
rections to the Maxwell dynamics that appear in a per- 
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turbative expansion of the effective action in QED in 2+1 
dimensions [37| ■ We focus on the extended Chern-Simons 
term which is Lorentz invariant, parity- violating, but no 
longer topological as the original Chern-Simons term. 



Extended Maxwell-Chern-Simons model 



which follows from Eq. (|44|) of the next section. 

By using the Eqs. ^ and (J7J) we can write down the 
conjugate momenta 



(13) 



Let us consider the extended MCS Lagrangian in 2 + 1 
dimensions in the Lorentz gauge 



£ = -\f^F" v + ^(DA a )(d p A,) - \(d,A»)\ 



(5) 



where g is a small coupling parameter. 

A direct calculation from the Lagrangian (JSJ) yields 

0^ = -^ + 9 r X ^~^An, (6) 



(7) 



The usual variation with respect to A a gives the gener- 
alized Euler-Lagrange equation 



d\ -^ A \ + o^dx- 



dA 



'd(8 x A 



1 d{d l _ l d x A a 



0, (8) 



and substituting Eqs. ([6]), ([7]) leads to the modified 
Maxwell equations 

d\F x<T + d"{d ■ A) + ^e a ^~'DF^ — 0, (9) 



which can be rewritten as 



e al3a d p ) UA a = 0. 



(10) 



Using the plane wave ansatz A„(x) — e M (fc)e , we 
find the two solutions 



fc 2 = 0, k 2 = \. 

9 



(11) 



Therefore, we see that the gauge field excitations are de- 
scribed by a massless and a massive mode (38[. Here we 
will use the convention rf- v ~ diag(l, — 1, — 1) together 
with the notation A; M = (fc°,k). 



B. The Hamiltonian formulation 

The Hamiltonian formulation for systems with higher 
time-derivatives was developed long time ago by Ostro- 
gradski [3(|. Accordingly, we consider A^(x) and A^x) 
as two independent configuration field variables with 
their corresponding conjugate momenta given by 



dA L 



Of 



dC 



(12) 



which read in components 

/" (d-A) I 



P l = F Oi + ^nA j + ^d F 0j) 



H° = fe^, 



We are using the conventions 

,012 _ 12 _ , 

e — e — + i, £oi2 — £12 — + L 



(14) 



(15) 



(16) 



(17) 



pOij _ ij 



1,3 = 1,2- (18) 



Now, considering the generalized Legendre transforma- 
tion Tic = P^A^ + li^A^ — C the canonical Hamiltonian 
density 7ic turns out to be 



He = P°A + -e^Pi Ilj - Pi(diA ) + -Fl - 
9 ± 9 

-n,(V 2 A0 + ie^^AoidiAj) + i(a ■ A) 2 , 

(19) 

where we have used P l Ai = 2e tJ PiUj / g — PidiA and 
Fl = 4H 2 /.g 2 . 

Examining Eqs. (Ti"4|) , (fi"6]) and (fTF)) we deduce the four 
primary constraints 



Xo 



n° - 9 - e^diAj , xi = P" + (d ■ A) + -.'u),X 



(20) 



Using the canonical Poisson brackets for the extended 
phase space variables 

{^(t,x),P I/ (t,x')} = r^ 3 (x-x'), 
{A lx (t,x),n v (t,x')} = r^ 3 (x~x'), (21) 

and after a straightforward calculation we find the non 
trivial elements of the algebra of constraints 

{Xi(x),Xo(x')} = 5 3 (x-x'), 

{^(x),xi(x')} = - 5 e^<5 3 (x-x'), (22) 

{¥>*(*), pV)} = 3 e^ 3 (x-x'). 
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Therefore, according to Dirac classification, we see that 
all the primary constraints fall into second class con- 
straints category. 

In order to search for additional constraints in the 
theory called secondary constraints, it is necessary to 
evolve the primary constraints with the total Hamilto- 
nian, which is defined adjoining the primary constraints 
X = {xo,Xi,<P 1 } as follows 

H T = J d 3 x (H c + <ToXo + trixi + <*V) , (23) 

where the variables {(To,fi,a ! } play the role of Lagrange 
multipliers. However, given that all primary constraints 
satisfy the weak condition 



{ X ,H T }^0, 



(24) 



there are no generations of secondary constraints in the 
theory. In addition, the Lagrange multipliers can be de- 
termined by comparing the Hamilton equations using the 
total Hamiltonian (|23| with the original ones of Eq. (fT0| . 
After some algebra we find 



oo = A , 



a x =0, a 1 = -d F { 



Oi- 



(25) 



In this way the total Hamiltonian is shown to be consis- 
tent with the equations of motion (|2"5|) and the set of con- 
jugate momenta (|T¥|) - (IT7jl . As an aside, let us mention 
that second class constraints can be imposed strongly to 
zero in which case we shall require the use of the reduced 
Hamiltonian Hr and the Dirac brackets given in the ap- 
pendix. 



C. The negative-energy mode 

To begin, let us express the general solution of Eq. (p~0|) 
as a sum of a massless and massive gauge fields as follows 



A^x) =A li {x) +B ll {x), 



such that A^(x) and B^(x) satisfy the equations of mo- 
tion 



UA a = 0, 



(27) 
(28) 



Now, let us expand both gauge fields in plane wave solu- 
tions 



A„(x) = 



d 2 k 



'(27r)°2fc a=o 



][>( A >(k)( aA (k), 



-i(kQXQ — k-x) 



(29) 



B M (x) 



d 2 k 



(7j M (k)b(k) e -<(*o*o-k-x) 



(2tt) j k 
7y;(k)&*(k) e^otfo-k-x)^ 



(30) 



where fco = |k| and fco = Wl + <? 2 |k| correspond to the 
massless and massive positive frequencies, respectively, 
see As in the usual case we are free to choose the 
massless polarization vectors to satisfy any normal- 
ization condition. Hence, we take 



eW(k). E l A '->(k) =v 



(31) 

However, according to Eq. (|28|) the massive polarization 
vector is uniquely determined and has necessarily to sat- 
isfy the normalization condition 

77(k)-7f(k) = -1. (32) 

One finds for both polarization vectors the relations 

e« 77, (k) r/* (k) = -igko, 7 7 (k) • rj(k) = 0, (33) 
rj(k)-k = 0, »fc(k) = Tfc(-k). 

Let us write the canonical Hamiltonian density in the 
form 



1 



(34) 



where the linear part in g is given by 



5Hc 



9 

+ 9 -e^A tl V 2 dpA 1 . 



(35) 



After a lengthy manipulation of replacing the gauge fields 
in Q34[) and using the properties of the polarization vec- 
tors (j3"3"|) the mode decomposition of the canonical Hamil- 
tonian is 



( 26 ) H c = / d 2 k(fc [a* (1) a (1) + a* (2) a (2) - a* (0) a (0) ] - k Q b*b). 



(36) 

We observe the negative-energy contribution coming 
from the time component of the massless gauge field in 
the quantity 

a *(0) o (0) 

which may be removed for instance 
switching to the physical Coulomb gauge as performed in 
the usual Maxwell theory. However, irrespective of which 
gauge is chosen, the negative-energy contribution from 
the massive mode will persist making the total energy 
unbounded from below. At this level a revision of the 
degrees of freedom is straightforward, see [39| . 



III. A BACKGROUND REVIEW 

In this section, the emphasis will be both to recall the 
canonical formulation for higher derivative field theories 
and to introduce the perturbative method according to 
the Refs. 0, [3l|, [HJ that will be used later on. 
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A. Hamiltonian formulation for higher derivative 
field theories 

Let us start with the action 



d 4 x C ((f>, d^cj), d^dvj), ... , d lil d lia ...d u 



(37) 



which for simplicity we assume to be a function of a scalar 
field and to depend on a finite number of N derivatives. 
A general derivation allowing for more general types of 



fields and infinitely many derivatives can be found in [40( • 
The standard procedure of extremizing the action (|37|) 
leads to the generalized Euler-Lagrange equation 



8C 



8C 



dC 



+(-1)^ 



d{d tl d u 4>) 
8C 



d(d^ ...d M 



0, (38) 



and analogously to the generalized energy-momentum 
tensor 



-sac- 



dC 



d v <j) ■ 



8C 



d{d^d^<j)) 



8C 



8L 



> 8(8,^8, 



Band, 



8L 



■■■ + (-l) N - 1 

xd m _ 1 d u <f>+ ... +(-l) N ~ 1 



8C 



d " 1 "' d, " f - 1 'd(d ltl ...d m _ 1 dtf 



8C 



■••d liN _ 1 d lt 

Using the Euler-Lagrange equation together with 

/ dC \ „ „ , / 8C 



d„4> -d^... 8^ N _ 2 
8 8 8 



\8(8 ai (j)) 



8 ai 8 v 



8(8 ai d a2 (j)) 



dai da.2 



it can be proven that the energy-momentum tensor sat- 
isfies the conservation equation d^T^ = 0, for a detailed 
derivation see Ref. [icj . 

The corresponding canonical formulation is given 
by a phase space of dimension 2N per space 
point, been characterized by N configuration 
field variables which by introducing the notation 
can be written as Qo(t,x) = 4>(t,x),Qi(t,x) = 



J 



8L 



d(8 ai . ■■d aN _ 1 c 



8L 



(39) 



</.W(t,x),...,Q iv _ 1 (f,x)=0( JV - 1 )(t ) x) ) where 

c>^(t,x) 



8t n 



(40) 



together with their corresponding conjugate mo- 
menta Po(t, x), Pi(t, x), . . . , P/v_i(t, x), and the canoni- 
cal Hamiltonian density which is identified with the com- 
ponent Tq, namely 



He 



-C 



8C 



d(8 (b) 



8C 



8C 



9 O 



8C 



dC 



d^8o(, 



(41) 



8C 



■■■ + (-l) N - 1 
+ (-l) Ar - 1 



8C 



d(8 fll ...d f , N _ 1 d 



8(8^ . . .8 fJ , N ^ 1 8 (l)) 
d^i ■ ■■d l j, N _ 1 d , 



<V - • ■ d tlN _ 1 — - ——- d (j) -8^... d m _ 2 — -—- d m _ x do<l) + ... 



8L 



<9(<9 Ml • ■■d tlN _ 1 d c l 
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The strategy to determine the conjugate momenta in the 
above expression will be to perform successive spatial in- 
tegration by parts in the fields to subsequently regroup 
the quantities multiplying the pure time-derivative terms 
8q (f>. Therefore, after successive integration by parts, the 
canonical Hamiltonian density (|4"Tj) can be given in terms 
of space phase variables as 

He = P (x)Q (x)+Pi(x)Qi(x)+... 

+PN-i(x)Qjt-i(x) - £, (42) 

from we read the explicit expressions for conjugate mo- 
menta 



P N -i(t,x) 



dC 

dQ N ' 



(43) 



Pi(t,x) 



9Qi+i fri + 

dC \ dP i+1 

dd mi . . .d mj Qi+i) dt 



(44) 



the index i runs trough i = 0, . . . , N — 2 and d mk stands 
for spatial derivatives. We have the compact expressions 
for the canonical Hamiltonian and the generalized sym- 
plectic form 



n(t) 



/d 3 x (^J2 P iQi- £j > ( 45 ) 



dPi(x,t) AdQ,(x,i) . (46) 




The basic Poisson brackets which can be read off from 
the above symplectic form are 



{Qifax^Pjfax!)} = %<5 3 (x-x'), 
i,j = 0,...,N-1. 



(47) 



Let us remark that due to possible degeneracy of the 
higher derivative theory, there may not be a unique solu- 
tion expressing Qi (x, t) in terms of the other canonical 

variables, occurring when ( ) = 0. In this case 

one can follow the Hamiltonian treatment for singular 
higher derivative systems developed in Ref. [4lj . 



B. The perturbative method 

Let us consider a framework of a system having a La- 
grangian density of the form 



(48) 



where C\ is an arbitrary function of the fields (not neces- 
sarily quadratic) that contains all the higher order deriva- 
tive dependence. In the limit of the coupling parameter 



g —* we recover the ordinary derivative theory defined 
by Co- 

The perturbative method for a canonically formulated 
higher derivative theory can be described in the following 
steps: 

(i) In order to obtain the appropriate Hamiltonian to 
order g p one starts by iteratively solving the equation of 
motion to the order 

(ii) Using the above iterated equation we express all the 
time-derivatives of the fields <j)^ (t, x) with q > 2 in terms 
of the first-order variables <p(t, x), (j>(t, x), called from now 
velocity phase space variables, and their possible spatial 
derivatives. This will introduce further contributions in 
powers of the perturbation parameter g which need to be 
maintained only up to the required order. 

(iii) We substitute the above expressions 
(iz) (i,x), q > 2 into the conjugate momenta Pj, 
the canonical Hamiltonian and the generalized symplec- 
tic form f2(£). This allows to recast up to the order p 
in terms of velocity phase space variables the first-order 
expressions 



H (p) 



(p) 



P 



O), 



+ o(g p+1 ) 



d^n^H^ <j)) + o(g p+1 ), 



(49) 



(50) 



n {p) (t) = J J d 3 xd 3 yf2W(i,x,y) 

x (#(*,x) A d0(t,y)) +<D(gP +1 ), (51) 

From above we can deduce the Poisson bracket 

{^,x),#,y)} = (fiW(t,x,y))- 1 , (52) 

the hat in Cl^ pS> (t, x, y) denotes possible dependence on 
spatial derivatives. 

(iv) In order to diagonalize (|5T|) . we search 
for an invertible change of variables from velocity 
phase variables (<f>(t, x), <j)(t, x)) to new canonical ones 
(Q(4>,4>), P{<j),4>)) in such a way that the new Poisson 
bracket {Q(t, x) , P(t, y)} are canonical to the order con- 
sidered. That is to say 



{Q(i,x), P(t,y)} 



d3z fJSQ(t,x)SP(t,y) 
\S(/)(t,z) 6(/>(t,z') 

SQ(t,x) SP(t,y)\ 



6j>(t,z') 8<f>(t,x) J 

x{cf>(t,z),4>(t,z')} 
= 5 3 (x-y)+0(g p+1 ) 



(53) 



(v) Having the kinematical theory completed we pro- 
ceed with the dynamics of the theory. Thus, replacing 
the new variables we rewrite the Hamiltonian as 

ff<*> = J d 3 xH^(4>(Q,P),<i>(Q,P)), (54) 
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and using the Hamilton equation 



(55) 



we write conjugate momenta P in terms of Q and Q. 
Finally, by using the Legendre transformation we arrive 
to the effective Lagrangian density 



C = P{Q,Q)Q-H{Q,Q). 



(56) 



Therefore, in terms of the new variables Q, P both the 
Hamiltonian density H(Q,P), together with the Poisson 
bracket |<3(i, x), P(t, y)| = <5 3 (x — y) define the physi- 
cal approximation of the system to the order considered. 
The first-order Hamiltonian will be bounded from below 
provided the initial one obtained from Co is. One can 
check that the Euler-Lagrange equations reproduce those 
of the original system to the order considered; a proof of 
self-consistency to all orders is provided in a mechanical 
setting in Ref. [33]. 

IV. THE EFFECTIVE THEORY 

In this section we describe the main ingredients and 
results of the application of the perturbative method to 
the extended MCS theory [I^]. The method has been 
previously applied for fermions and scalars in Refs. [H, 

A. The reduced phase space 



Replacing into the set of canonical variables ([14]) we ar- 
rive to 



P l = F M + ~e ij V 2 Aj - 9 -e^d 3 A + 0(g 2 ), (59) 



p° = -(s-A)-| c «a i ^ > 



H° = ZeVdiAj, 



n' 



L e Um 



0.; • 



(60) 
(61) 
(62) 



The perturbed Hamiltonian density l|3~4")) gets expressed 
by 



n 



(i) 



■\{A 2 - A V 2 A ) + ±(A 2 AiV'Ai) 



/y ( ''<V,.i,,i :: • /y.'-'V-M,,!, + !l ('u),.\X- ! .\, ! 

+0(g 2 ), (63) 
and the perturbed symplectic form 



n (1) (t) = J d 2 x (dAi A dA z - dA Q A dA 

-g€ lj [d{d l A 3 ) A dA + didjAo) A dA t ] 

+ ^e ij [dAi A dAj - d{V 2 Ai) A dAj]) + 0{g 2 ). 

(64) 



We start by iteratively solving the equation of motion 
(fT0| to the lowest order approximation (p = 1). A first 
iteration gives 



UA a = g\d a U{d ■ A) - D 2 A a ) « 0(g z ) 



(57) 



Therefore, to the order considered, our iterated equation 
of motion is the basic equation 



Let us consider the equivalent two-point split expression 



n( 1 )(t) = i J J rf 2 x d 2 x'Q^(t;x,x') dz a {t,x) Adz b (t,yi'), 

(65) 



A a = V 2 A a +0(g 2 ). 



(58) 



where 





" 








1 


-gd 2 


gdi ' 










-<?v 2 


gd 2 


-1 





^;x,x') = 




-1 


gd 2 




-gdi 


-gdi 








-l 






-gd 2 


l 











g 




. gdi 





1 





-g 






<S 3 (x-x'), 



(66) 



and the notation is such that z a — (A , A 1 , A 2 , A , A 1 , A 2 ) with a — 1,2,3,4,5,6. Notice that a minus sign has 
appeared in the linear derivatives due to an integration by parts. Also, here and in the following, to avoid ambiguities 
we will consider the action of the derivatives with respect to the unprimed coordinate x. 
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The inverse symplectic matrix is easily computed to give 



(ft£(t;x,x'))- 



" 








-l 


gd 2 


-gdi ' 










9 


-gd 2 


l 









1 


-g 

-gd 2 




gdi 


gdi 







l 




<5 3 (x 


gd 2 


-1 











-gV 2 




. -gdi 





-i 





gV 2 








x'). 



(67) 



It can be checked that the usual properties for the ma- 
trices (|66l) and (|67p . such as the antisymmetry property 



^ 1 h ) (t;x,x') = -^ ) (t;x',x) 
together with the inverse relation 



(68) 



d 2 x' n£(t; x, x')^ (t; x', x"))- 1 = <5 ac <5 3 (x - x"), 



are satisfied up to linear order in g. 
Now, considering that 



(69) 



{**(t,x), z 6 (t,xO} = (oi 1 6 ) (t;x ) x'))- 1 , (70) 
from (|67p the nonzero brackets are 

{A)(*,x),A (t,xO} = -<5 3 (x-x'), 
{A^x), x')} = c%<5 3 (x-x'), 
{^(t,x),A,-(t,xO} = 5 3 (x-x'), 

{A(«,x),i (t,x')} = -^^-^(x-x'), 

{i4o(*,x), i^x')} - ^^^(x-x'), 
{i^,x),A,(i,x')} = -. g e^V 2 <5 3 (x-x'). (71) 



B. New canonical variables 



or alternatively 

(A(t, x), i(t, x)) , 5r„ (,4(t, y), y)) } 

= ^<5 3 (x-y)+0(. g 2 ), (75) 

which may be verified with the use of (fTTj) . It is important 
to mention that the above realization can be achieved by 
more than one set of canonical variables. Hence, in order 
to fix a unique set of canonical variables we will impose 
additional criteria on the gauge field redefinition: These 
requirements are: (i) we choose the new gauge field 
to transform covariantly under Lorentz transformations 
and (ii) to satisfy the same gauge fixing condition as the 
old one. 

Here the simplest way to proceed is to select a new 
set of canonical variables and check if Eqs. (f?4"|) , ([75)1 to- 
gether with the two above conditions are verified. There- 
fore, let us consider the twisted gauge redefinition 



which in components becomes 



A = Ao + ^diAj. 



(76) 

(77) 
(78) 



We define the corresponding canonical momenta by 

n = -(d-A), (79) 



Until now we have the effective theory given in terms 
of velocity phase space variables with their corresponding 
Poisson brackets exhibiting the non standard form (|7Tj) . 
According to the method we need to search for a new set 
of canonical variables for which the Poisson brackets (|7"Tj) 
are diagonal. The effective theory will be constructed in 
terms of these variables. Therefore, let us consider the 
invertible change of variables 

A - A{A,A), (72) 

A -» w(A,A), (73) 
such that (|64|) can be put in the form 

fiW(t) = [ d 2 xdn f *(t,x)AdA ll (t,x)+0(g 2 ), (74) 



5P = F 0i + |e«V% - |e^i - ^didiAm. (80) 

One can check that the new gauge field respects the gauge 
fixing condition 

XG .F = (d-A)(d-A), (81) 

and transforms covariantly under Lorentz transforma- 
tions, which we further prove. 

It is not a difficult task to check that the new vari- 
ables diagonalize the Poisson brackets, using (|7Tj) let us 
compute the Poisson brackets 

{2 (*,x),5fo(i,xO} = {a + ^diAj,-(d-A)} 

= -{A c + ^ i a i A j ,Ao-d k A k } 

^ 3 (x-x')+dV), (82) 
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{^(i,x),F(t,x')} = [Ai + £ e im F 0m ,F 0j + 9 -e^ k V 2 A k 



9 -e^ m d m A Q - |e iro c^ m } 



(83) 



and 



{Mt,x), Mt,x')} = [A + 9 -e im F Qm ,A j + f e^Jbm} 

5 



9e y + ^ m S im - Y m 8 mj + 0{g 2 ) 

(84) 



= 0, 

from where we finally arrive to 

x),Sv(t,x')} = r^,5 3 (x-x')+0(.g 2 ). (85) 

The same result follows by considering the inverse change 
of variables 



A = Ao-^d.A,, 



Ai — Ai 



-€ J TT J , 



A Q = d k A k -n -^diP, 

Ai = tt* + d t A + ^' *<>,<>, A, - |e^o - |e y V 2 !,, 

(86) 



and replacing in the perturbed symplectic form 

n«(t) = J d 2 x (dAi A dA, - di A dA 

-ge^ididiAj) A (L4 + d(d j A Q ) A dA^ 
+ ^[dAi A dij - d(V 2 ^) A dAj]), 

which after some algebra produces 

n^{t) = J d 2 x fi^t.xjAdA^i.x), 
in agreement with the result 

C. The low-energy effective Lagrangian 



Recall from the Eq. (|6"3"|) the canonical Hamiltonian 
density written in the old variables 

= -^ + ^ V 2 ^o + ii 2 -^V 2 ^ 

-ge^'^ijio + ge ij V 2 AjAi + ge^ d l A j V 2 A 0l 

(89) 

which with the use of the inverse transformations has 
to be expressed in terms of the new canonical variables. 



Using the inverse transformations (|86|) we have the ex- 
pressions 

-\Al = -\{d k A k f l -{H f +MdkA k ) 

-i(a^o) 2 = -\{d % Av) 2 + ^(didiA^idiAo), 
±A 2 = l(^) 2 + l(d t A ) 2 +n*(d t A ) 

+ ^{^d 3 d k A k - |"/U„ - fe^V 2 !,-) 
- 1^4, V 2 !,-, 



and the ones linear in g 



(9 i ^) 2 = ^(^ i ) 2 -| e ^^l ; . 



(90) 



.7<''/U,.l, = -ge y d t F(cM fe )+ge l W J )^>, 
//<''V-'.l.,,i, = -y/'-' V J .l, • //( '-'V J .UrU.,;. 
.^'•'^.TV-.l, = .ge^(cU,)V 2 lo. (91) 
Replacing in the Hamiltonian, we arrive to 



= -^ 2 + ^(^) 2 +^o(<9 fc A fc )-A (a 4 P) 
1 



(92) 



where Fij(A) = 9j^4j — <9jAi. 

Considering the Hamilton equation A^ = j^^ii^ 1 ) j 

we obtain the canonical momenta in terms of configura- 
tion variables and their first time derivatives 



n = -(d-A), = Ai — diA Q . 



(93) 



The effective Lagrangian is obtained via the Legendre 
transformation 

C = ^ l A t ,-H {1) (A,n). (94) 
Finally, using the relations we obtain 

C = -- A F, v (A)F^(A)- l -{d-Af. (95) 

Therefore, considering that A^ — A^ + |e MQj a d a A l3 and 

F liV {A)d ll A v — dyAy, the Maxwell theory is correctly re- 
produced by setting g — > 0. 

The standard variation produces the equation of mo- 
tion 



UA a = 0, 



(96) 



and by considering plane waves with respect to the new 
gauge field gives the solution k 2 = 0. We see that only 
the massless solution has been retained in the effective 
field theory while the massive mode has been removed. 
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1. Gauge invariance 

Recall that classically appropriate conditions can be 
imposed on the boundary data in order to set d ■ A = 
everywhere so we can lift this condition from the effective 
theory. Therefore, let us consider the Lagrangian 



In order to find expressions analogous to those of Eqs. 



C = -\f i1v {A)F^{A), 



(97) 



which is manifestly invariant under the gauge transfor- 



mation A'^ — A^ — d^A. 



2. Lorentz symmetry 

Lorentz invariance follows naturally in the effective 
theory by proving covariance of the gauge field A 11 . 
Therefore, consider the Lorentz transformation 



A" -> A" = A" + ^ vp d v ,A p >, 



(98) 



and replace the transformation for the original fields 



(99) 



together with the relation e" Vp ' A a v , A } = e aaX A^. This 
gives the desired relation 



A" = At 



9 -e a ^d a A x ) = A',', .1". (1(H)) 



by 



CVT symmetries 



Parity transformation V is defined in 2 + 1 dimensions 



x° - x\ 



^2 ^ x 2 



(101) 



which corresponds to a reflection in just one of the spa- 
tial axes [44| . Notice that this transformation yields the 
improper transformation defined to have det A = — 1, in- 
stead of the space inversion x — * — x one is familiar in 
three spatial dimensions. 

The original gauge field transforms as 



VA°{x°,x)V~ 1 = A°(x°,x'), 

VA 1 {aP,x)V- 1 = -A 1 (x°,x')> 
VA 2 (x°,x)V- 1 = A 2 (x°,x'), 

where x = [x , x 2 ) and x' = (—a; 1 , x 2 ). 
From (|76p we have in components 



(102) 



A°(x°,x) = A°(x°,x) + I (diA 2 (x°,x) - d 2 A 1 (x°,x)) , 
A^x ,*) = A 1 (x°,x) + | (a 2 A (x°,x) - d A 2 (x°,x)) , 
A 2 (x°, x ) = A 2 (x°,x) + | (flbAi(ar ,x) - ^(x^x)) . 



(|102[) it is convenient to set A^(x) = e M (fc)e 
Eq. (fT03l) becomes 



A»(x°,x) = T»(k°,k)A a (x Q ,x), 



such that 



T(fc°,k) 



1 ~~igk 2 /2 igki/2 
—igk 2 /2 1 —igko/2 
igki/2 igko/2 1 



Under a parity transformation we find 

VA^(x°,x)V- 1 = R p a (k Q ,k)A a (x ,x') 

with 



R(k°,k) 



1 igk 2 /2 —igki/2 
—igk 2 /2 —1 —igko/2 
—igki/2 —igko/2 1 



where 



(104) 



(105) 



(106) 



(107) 



We need the original gauge field evaluated at x', so from 
(fTUI)) we get 

A"(a: ,x') = T p a (k ,k)A a (x°,x'). (108) 
Replacing (fTU8"|) in (fTUH|) we obtain 

VA"{x°,x)V- 1 = T^(fc°,k)i Q (x°,x'), (109) 

where the matrix of the transformation T(k , k) = 
Rik^tyT-^k ^) is find to be 



1 igk 2 -igh 
-igk 2 -1 -igk 
-igki -igk a 1 



i -g& 2 gd[ 

gd' 2 -1 gd'o 
gd[ gd'o 1 



(110) 

From (f9"5"|) and (| 1 10[) we have that the Lagrangian trans- 
forms as 



VCix^- 1 =C(x') + 0(g 2 ) 1 



(111) 



(103) 



where we have used T T r/T = rj + 0(g 2 ). We conclude 
that the effective theory is invariant under parity trans- 
formations to the order considered. One can in addition 
check that the effective Lagrangian is time reversal T and 
charge conjugation C invariant. 



V. CONCLUSIONS 

In this paper we have studied higher derivative field 
theories within a framework in which their corresponding 
higher derivative operators are regarded as small correc- 
tions to an ordinary-derivative field theory. Historically, 
the use of higher derivative operators have been avoided 
in field theories due to the many problems encountered 
with their formulation. Perhaps the most significant ones 
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are the proliferation of extra degrees of freedom with re- 
spect to the ordinary ones therefore going against the 
premise that corrections should introduce small devia- 
tions together with the appearance of Hamiltonians be- 
ing unbounded from below. Indeed, these problems arise 
regardless of whether these corrections appear as free or 
interacting terms in the Lagrangian. In general the so- 
lutions in the higher derivative theory can be classified 
according to their analytical behavior in the limit — >• 0. 
We associate the analytical solutions to low-energy de- 
grees of freedom and the nonanalytical solutions to high- 
energy degrees of freedom. In this regard, we see that 
the appropriate way of treating such corrections as small 
perturbations is to implement techniques in order to con- 
trol and subsequently remove the high-energ y d egrees of 
freedom. The perturbative formulation [{J, l3l| accom- 
plishes this by suppressing the excitation of high-energy 
modes in a way consistent with the exact evolution and 
only allowing further excitations of the low-energy modes 
already present in the zeroth order theory. 

The main goal in this work has been to test the symme- 
tries in the low-energy regime of higher derivatives field 
theories. In particular we have focused on the higher 
derivative Chern-Simons theory. It is well known that 
symmetries in an effective theory will depend on the 
scheme of approximation used to obtain the low-energy 
limit. For example, the removal of all nonanalytical terms 
in the Hamiltonian (|19p justified within an analytical ex- 
pansion would lead to the violations of Lorentz invari- 
ance [35|. Also, we would have parity invariance to all 
orders with no g dependence by applying the equations 
of motion directly on the extended MCS Lagrangian (36| , 
which follows by replacing the lowest order iterated so- 
lution OA a rs in the original Lagrangian. 

To summarize, we have developed the Hamiltonian 
formulation for the higher derivative MCS theory in 
2 + 1 dimensions. Given the singular nature of the sys- 
tem we have been required to follow the Dirac method 
to impose second class constraints strongly from where 
we have computed both the Dirac brackets and the re- 
duced Hamiltonian. In addition, we have exhibited in 
the Hamiltonian the negative-energies producing insta- 
bilities in the interacting case. The perturbative method 
was implemented for the higher derivative Chern-Simons 
Lagrangian and an effective field theory that describes 
corrections to the low-energy physics with a Hamilto- 
nian well bounded from below and low-energy solutions 
were constructed. We have introduced a prescription for 
the choice of the new canonical field variables which are 
unique according to the following criteria: (i) the new 
gauge field transform covariantly under Lorentz trans- 
formations and (ii) it satisfies the same gauge fixing 
condition as the old gauge field. As suggested above, 
prior to (i) and (ii) there is some arbitrariness in the 
choice of new canonical variables since the diagonaliza- 
tion of the symplectic form may be achieved by differ- 
ent sets of canonical variables. Nevertheless, since they 
are connected by canonical transformations they all de- 



scribe physically equivalent theories. Although it is be- 
yond the scope of this paper, by adopting similar argu- 
ments than before we can provide an expression for the 
new gauge field to the next order which is uniquely given 
by = + Bg& a Pd a Afi + ng 2 UA^ with k, some 
real numbers. Finally, we have obtained an effective field 
theory which is parity invariant to the order considered 
therefore recovering symmetries one could expect at low 
energies. 
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APPENDIX A: DIRAC BRACKETS 

In this appendix we follow Dirac procedure to re- 
duce second class constraints from the higher deriva- 
tive theory. To begin, let us introduce the notation 
Pa(xoi Xii <f l )i with A = 0,1,1,2. The matrix of the 
second class constraints will be denoted by 

CU B («;x,x') = {p A (i,x),p JJ (t,x')}. (Al) 

From (f2"2"| we have 





- 


-1 





" 


Cab = 


1 





-gd 2 


-gd 2 



gdi 

9 




.0 


gdi 


~9 






£ 3 (x-x'). (A2) 



The inverse matrix is 





- 


1 


-dx 


-da 1 


^AB ~ 


-1 

-di 












-1/5 




_-d 2 





1/5 






The nonzero components are 
C 5I 1 (x,x') = 5 3 (x-x'), 

C 5 - i 1 (x,x') = -9^ 3 (x-x'), (A4) 

C i - 1 (x,x') = --e 4 ^ 3 (x-x'), i,i = l,2. 
5 

And as usual Dirac brackets are defined as 

{X, Y} D = {X, Y} - {X, y A } C A X B { VB , Y}. (A5) 
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After some calculation the nonzero Dirac brackets are and the reduced Hamiltonian density is 

{A (t,x),io(t ) x')}D = -5(x-x , ) ) . 2 / .. EL a o 

{Mt,x),P (ty)} D = S(x-x'), K R =P A - -m (e« Pj + -i + §V 2 A 

{4.(t,x),p*(ty)}2, = |e«a, «(x - x'), _ g iAo ( fl + £ e «V 2 ^) + 7*8 + A) 2 . 

V 2 / 4 J 2 

{i t -(t,x) J i i (t,x')}j 3 = --e^(x-x'), (A7) 

{Ai(t, x), Ao(t, x )}■£) = — 9j5(x— x ), Finally, the reader may check that the reduced Hamil- 

t i , . _ . ,., l^c/ / \ tonian together with the Dirac brackets give the exact 

{A i (t,x),P (i,x')}.D = -didhc-x'), , . , , . ,. , 6 

/J 2 higher derivative equations oi motion. 

{i i (t,x),P^(t,x / )} 1J = I^S^x-x'), 
{i 4 (i,x),n^(t,x / )} 1 , - i<%«5(x-x'), (A6) 

{i (t,x),P l (i,x')}D =9 i <5(x-x / ), 
{P (i,x),ff(i,x')} D |6 y 9^(x-x'), 

{n (t,x),P ! (i,x')}D = -fe^(x-x'), 

{P i (i,x),P^"(t,x')}D|e iJ V 2 ( 5(x-x'), 

{n*(i,x),rP(i,x')} D - -|e^(x-x'), 
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